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We give a general method for finding an exact analytical solution for the two state curve crossing 
problem. The solution requires the knowledge of the Green's function for the motion on the un- 
coupled potential. We use the method to find the solution of the problem in the case of parabolic 
potentials coupled by Gaussian interaction. Our method is applied to this model system to calcu- 
late the effect of curve crossing on electronic absorption spectrum and resonance Raman excitation 
profile. 
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I. INTRODUCTION 



Nonadiabatic transition due to potential curve (or sur- 
face) crossing is one of the most important mechanisms 
to effectively induce electronic transitions in collisions 
[T^ |. Two state curve crossing can be classified into the 
following two cases according to the crossing scheme: (1) 
Landau-Zener (L.Z) case, in which the two diabatic po- 
tential curves have the same signs for the slopes and (2) 
non-adiabatic tunnelling (N.T) case, in which the dia- 
batic curves have the opposite sign for slopes. There 
is also a non-crossing non-adiabatic transition called the 
Rosen-Zener-Demkov type P, 01; in which two adiabatic 
potentials are in near resonance at large R. The theory 
of non-adiabatic transitions dates back to 1932, when the 
works for curve-crossin g an d non-crossing were published 
by Landau [l^, Zener |l6l | and Stueckelberg and by 
Rosen and Zener [ll] respectively. Two categories might 
be classified for finding an exact analytical solution of the 
curve crossing problem. The first is that an exact analyt- 
ical solution can be obtained for the whole region of the 
variable (say x here, see in the next section). For exam- 
ple, Osherov and Voronin solved the case where two di- 
abatic potentials are constant with exponential coupling 
(igj . C. Zhu solved the case where two diabatic poten- 
tials are exponential with exponential coupling [20j. In 
our earlier publications we have reported the analytical 
solution in the case, where two or more arbitrary poten- 
tials are coupled by Dirac Delta interactions [2l|-|26|. The 
second is that an exact analytical solution is only possi- 
ble for the asymptotic region. Then, physical quantities 
such as eigenvalues, scattering matrices can still be solved 
in an exact analytical form, providing that the connec- 
tion problem of the asymptotic solution is known. The 
Stokes phenomenon [27] of asymptotic solution of the or- 
dinary differential equation provides a. powerful tool to 
deal with these kinds of problems [28l - t30 |. Generalizing 
the real variable to the complex variable and tracing the 
asymptotic solution around the complex plane, the con- 
nection matrix which connects the asymptotic solution 
in the complex plane can be expressed in terms of Stokes 
constants. Recent work by Zhu and Nakamura ^30,] found 
an exact analytical solution of the Stokes constants for 
the second-order ordinary differential equation with the 



coefScient function as the fourth-order polynomial. In 
this way, exact analytical solutions of scattering matri- 
ces were obtained for the two state linear curve crossing 
problem with constant coupling (sij . 



II. THE MODEL 

We consider two diabatic curves, crossing each other. 
There is a coupling between the two curves, which causes 
transitions from one curve to another. This transition 
would occur approximately in the vicinity of the crossing 
point. In particular, it will occur in a narrow range of x, 
given by 



\Vi{x)-V2ix)\ ^ \V{Xc)\ 



(1) 



where x denotes the nuclear coordinate and Xc is the 
crossing point. Vi and V2 are the diabatic curves and 
V represent the coupling between them. Therefore it is 
interesting to analyze a model, where coupling is local- 



ized in space near Xj given b> 
kj is the coupling strength [21|, 



■Xj), where 
The majority of 
the problems of interest, however do not corresponds to 
a localized coupling and one requires different forms of 
coupling V{x) for proper description of dynamics in dif- 
ferent cases. We express the arbitrary coupling function 
V{x) in terms of a linear combination of Dirac Delta func- 
tions [H, [111 ■ Expressing the arbitrary coupling function 
V{x) in terms of Dirac Delta functions has the advantage 
that it can be solved exactly by using analytical methods 
[2l| - [2^ . An arbitrary coupling V{x) can be written as 



Vix) 



dx'Vix')6{x ~ x') 



and the above integral can be discritized as 



N 



= '^kjS{x-Xj), 



here kj are constants, given by 



WjVixj). 



(2) 



(3) 



(4) 



The weight factor Wj varies depending on the scheme of 
discritization used 1321 . 
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III. EXACT ANALYTICAL SOLUTION 

We start with a particle moving on any of the two dia- 
batic curves. The problem is to calculate the probability 
that the particle will still be on any one of the diabatic 
curves after a time t. We write the probability amplitude 
as 



ipi{x,t) 



(5) 



where ipii^x, t) and '>p2{x, t) are the probability amplitude 
for the two states, ^'(a;, t) obey the time dependent Schrii 
dinger equation (we take h = \ here and in subsequent 
calculations) 



where H is defined by 



H = 



Hi{x) V{x) 
V{x) H2{x) I ' 



where Hi{x) is 



1 92 

H^{x)=- — J-^ + V^{x) 

2m ox'^ 



(6) 



(7) 



(8) 



We find it convenient to define the half Fourier Transform 
'^{E) of *(<) by 



Jo 

Half Fourier transformation of Eq. ([6]) leads to 



(9) 



1,{E) \ _J E-H, -V 
ME)) [ -V E-H2 



MO) 
MO) 



This may be written as 

*(£;) = iG(£;)*(0). 



(11) 



(10) 



G{E) is defined by {E - H) G{E) = I. In the position 
representation, the above equation may be written as 

/oo 
G{x,xo;E)'^{xQ,E)dxo, (12) 
-00 

where G{x,xo; E) is 

G{x,xo;E) = {x\{E-H)-^\xo). (13) 

Writing 

G(x,xo;i?)=f^";"'"°'^^| G,,(.,.o;i?n (14) 
\G2i{x,xo;E) G22ix,xo;E) J 

and using the partitioning technique [3(| we can write 

G^^ix,xo;E) = {x\[E-Hi-ViE-H2)-'vM\xo). (15) 
The above equation is true for any general V . This ex- 
pression simplify considerably if V is expressed as a sum 
of delta functions [H, [s^ . In that case V may be written 
^sV = Ef=i K,S, = Ef=i K,\x,){x,\. Then 

TV 

G,,{x,xo;E) = {x\[E-Hi-J2K^GUx,,x,;E)Sj]-'\xo), 

(16) 

where 

G°{x,xo;E) = {x\{E~H2)-'\xo), (17) 

and corresponds to propagation of the particle starting 
at xq on the second diabatic curve, in the absence of 
coupling to the first diabatic curve. Now we use the 
operator identity 0, HH] 



N N N 

{E-H^-Y,Kpl{x„x,-E)S,)-^ = {E-H,)-^+{E-H{)-^Y.^p2{xj.x,-E)S,[E-H,-Y,Kpl{xj,xf,E)S,]- 

, 08) 

I 

Inserting the resolution of identity / = dy \y) {y\ in we obtain a set of linear equations, which can be written 
the second term of the above equation, we arrive at as 

N 

Gi,{x,xo;E) ^Gl{x,XQ-E) + Y,Kpl{x,Xj-E) (19) 
X G2{xj,Xj;E)Gii{xj,XQ\E). 



Considering the above equation at the discrete points Xi , 



(20) 
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where the elements of the matrices A = [a.y ], P = [pi] 
and Q = [qi] are given by 

a^j = -kfGl {xi , Xj ; £;)G'§ {xj , xj ; E) + d^j (21) 
Pi = Gii{xi,xo;E) 
qi = G'lixi^XQ-^E) 

One can solve the matrix equation i.e. Eq. (20) eas- 
ily and obtain Gii{xi,Xo;E) for all Xi. Eq. (19) then 
yield Gii{x,xo; E). Similar one can derive expressions 
for Gi2{x,xo; E), G22{x,xo; E) and G2i{x,xo; E). Using 
these expressions for the Green's function in Eq. (11) we 
can calculate '^{E) explicitely. The expressions that we 
have obtained for '^'{E) are quite general and are valid 
for any Vi{x) and V2{x). 

IV. ELECTRONIC ABSORPTION SPECTRA 
AND RESONANCE RAMAN EXCITATION 
PROFILE 

In this section we apply the method to the problem 
involving harmonic potentials. We consider a system of 
three potential energy curves, ground electronic state and 
two 'crossing' excited electronic states (electronic transi- 
tion to one of them is assumed to be dipole forbidden and 
while it is allowed to the other) [s^, H^- We calculate 



In the above equation Hmb,ei{x) and Hmb,e2{x) describes 
the vibrational motion of the system in the first elec- 
tronic excited state (allowed) and second electronic ex- 
cited state (forbidden) respectively 

19^ 1 

i?m6,ei(a;) = - — ^ + -m£;l(a;-a)2 (23) 

and 

H,^b,e2{x)^-^g^ + ^rnEUx^bf. (24) 

In the above m is the oscillator's mass, Ea and Ep are 
the vibrational frequencies on the allowed and forbidden 
states and x is the vibrational coordinate. Shifts of the vi- 
brational coordinate minimum upon excitation are given 
by a and 5, and V12 (V21) represent coupling between the 
two harmonic potentials which is taken to be 

V21 (x) = 1^12 (x) = i^oe-"("-"=)' , (25) 
where Kq represent the strength of the coupling. The 



the effect of 'crossing' on electronic absorption spectra 
and on resonance Raman excitation profile. The propa- 
gating wave functions on the excited state potential en- 
ergy curves are given by solution of the time dependent 
Schrodinger equation 




Configuraiional Coordinate 



FIG. 1: Schematic diabatic potential energy curves that illus- 
trate the model. The forbidden state is labeled F; the allowed 
state is labeled A. 



(22) 

I 

intensity of electronic absorption spectra is given by (ssl . 

Ia{E) cx Re[J^^ dx dx,^T'* (^) 

iG{x,xo-E + iT)^T\xq)1 (26) 

where 

G[x,XQ;E+iT) - {x\[{Eo/2+E-Eeg)+iT-H,,b.e]-^\xo). 

(27) 

and 

_ / Hy,b,el{x) Kq\Xc){Xc\ \ ,„„s 
[Ko\x,){x,\ H,,,.,2{x) ) ^^^> 

Here, T is a phenomenological damping constant which 
account for the life time effects. ^"*^(a;, 0) is given by 

^T\^.^)=[^'t^)^ (29) 

where Xi{x) is the ground vibrational state of the ground 
electronic state, Eq is the vibrational frequency on the 



d f ri'\x,t) 



dt\r2'\x,t) 



Hyib.el{x) 
V2l{x) 



Vi2{x) 
Hvib,e2 (x) 



ri^'ix^t) 
r2'\x,t) 
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FIG. 2: Calculated electronic absorption spectra with cou- 
pling (solid line) and without coupling (dashed line). Here the 
values for the simulations are Eq = Ea = Ef ~ 400 cm~^ , 
r = 450 cm"\ £A = 10700 cm"\ ef = 11500 cm"\ 
Ko = 1200 cm"\ a = 0.2594 A"^ and Xe = -0.1991 A. 

ground electronic state, ea is the energy difference be- 
tween the excited (allowed) and ground electronic state, 
and for the forbidden electronic state it's value is sp- 
Similarly resonance Raman scattering intensity can be 
expressed in terms of Green's function and is given by 

iG{x,xo;E + iT)^';'\xo,0)\^. (30) 
In the above ^'J*^(a;, 0) is given by 

^^(^'O)^ (%^''^) ' (31) 

where Xfi^) is the final vibrational state of the ground 
electronic state. As G'f{x,Xo;E) for the harmonic po- 
tential is known ji^], we can calculate G{x,xo; E). We 
use Eq. (pOl) to calculate the effect of curve crossing on 
resonance Raman excitation profile. 

A. Results using the model 

In the following we give results for the effect of curve 
crossing on electronic absorption spectrum and resonance 
Raman excitation profile in the case where one dipole 
allowed electronic state crosses with a dipole forbidden 
electronic state as in Fig. [T] As in [3g], the ground 
state curve is taken to be a harmonic potential energy 
curve with its minimum at zero. The curve is constructed 
to be representative of the potential energy along a 
metal-ligand stretching coordinate. We take the mass as 
35.4 amw and the vibrational wavenumber as 400 cm~^ 
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FIG. 3: Calculated resonance Raman excitation profile for ex- 
citation from the ground vibrational state to the first excited 
vibrational state, with coupling (solid line) and without cou- 
pling (dashed line). Here the values for the simulations are 
Eo = Ea = Ef ^ 400cm-\ F = 450cm-\ £a = 10700cm-\ 
EF = 11500 cm~\ Ko = 1200 cm-\ a = 0.2594 1"^ and 
= -0.1991 A. 



[38| for the ground state. The first diabatic excited state 
potential energy curve is displaced by 0.2 A and is taken 
to have a vibrational wavenumber of 400 cm~^. Transi- 
tion to this state is allowed. The minimum of the po- 
tential energy curve is taken to be above 10700 cto~^ 
of that of the ground state curve. The second diabatic 
excited state potential energy curve is taken to be an 
un-displaced excited state. On that potential energy 
curve, the vibration is taken to have same wavenumber of 
400 cm^^. Its minimum is 11500 cm^^ above that of the 
ground state curve. Transition to this state is assumed 
to be dipole forbidden. The two diabatic curves cross at 
an energy of 11765 cm~^ with Xc — —0.1991 A. Value of 
Kq we use in our calculation is Kq — 1200 cm^^ and the 
value of a we use in our calculation is a = 0.2594 A~^. 
The lifetime of both the excited states are taken to be 
450 CTO~^. The calculated electronic absorption spectra 
is shown in Fig. The profile shown by the dashed line 
is in the absence of any coupling to the second potential 
energy curve. The full line has the effect of coupling in 
it. The calculated resonance Raman excitation profile is 
shown in Fig. [3] The profile shown by the full line is 
calculated for the coupled potential energy curves. The 
profile shown by the dashed line is calculated for the un- 
coupled potential energy curves. It is seen that curve 
crossing effect can alter the absorption and Raman ex- 
citation profile significantly. However it is the Raman 
excitation profile that is more effected. 
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V. CONCLUSIONS 

We have proposed a general method for finding the ex- 
act analytical solution for the two state curve crossing 
problem. Our solution is quite general and is valid for 
any potentials for which Green's functions for the mo- 
tion in the absence of coupling is known. We use the 
method to find the solution of the problem in the case 
of parabolic potentials coupled by Gaussian interaction. 
Our method is used to calculate the effect of curve cross- 
ing on electronic absorption spectrum and on resonance 
Raman excitation profile. We find that Raman excita- 



tion profile is affected much more by the crossing, than 
the electronic absorption spectrum. 
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